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La pr6sente invention conceme une methode pour compacter et filtrer des donn^es 
sismiques dans le but d'en extraire au mieux des informations sur la nature du sous-sol. 

Cette m6thodologie s' applique k la compaction et au filtrage des jeux de donn6es 
sismiques "multicubes", i.e. soit des donn€es sismiques avant sommation issues d'une 
5 campagne sismique 3D regroup6es en un certain nombre de cubes de donnees 
conrespondant chacun par exemple ii un mSme d6port (cubes de traces sismiques dits iso 
offsets) ou k un mgme angle d'incidence (dits cubes iso angles d'incidence), ou i des cubes 
d'attributs sismiques obtenus par traitement des amplitudes sismiques comme, par 
exemple, des cubes de propri6t^s 61astiques issus d'une inversion stratigraphique jointe de 

10 donn6es sismiques avant sommation, ces cubes 6tant obtenus par la mise en oeuvie de la 
mdthode d6ciite par exemple dans le brevet ER 2 800 473 (EP 1 096 270) du demandeur. D 
s'agit alors de mettre en forme I'information sismique avant sommation contenue dans des 
portions de traces sismiques form6es par exploration d'une zone souterraine, k I'aide de 
. nouveaux attributs repr6sentatifs qui sont ensuite exploits pour I'interpretation geologique 

15 ou I'interpretation dHm reservoir du sous-sol. EUe peut aussi s'appHquer h. la compaction et 
au filtrage de donnees sismiques 3D avant ou apres sommation de campagnes repa^esV 
dans le temps (sismique 4D), ceci afin d'interprdter les variations de contenu en fluides^ 
d'un reservoir p6troIier en cours de production. V 

Etat de la technique 



20 



Les mesuies sismiques sont classiquement utilises pour foumir des informations 
compl6mentaires, par rapport aux donn6es de puits. sur les variations des propri6t6s 
g6ologiques des formations du sous-sol. L*interpr6tation sismique est bas6e sur 1' analyse 
d'attributs afin d'extraire ces propri6t6s g^ologiques k partir de I'information sismique, et 
ce dans une zone d'int6r@t d'une formation souterraine (zone r6servoir par exemple). 

25 Une approche standard pour 1' analyse d'attributs implique I'utiUsation de m^thodes 

statistiques de reconnaissance de forme et d'estimation afin de calibrer les attributs 
sismiques pertinents avec I'information disponible obtenue au droit de diff^rents puits 
for&. 



Par exemple. par le brevet FR -A 2 768 818 (US 6 051 651 du demandeur). on 
connait une methode permettant de detecter les classes potentielles dans une population 
d'6v6nements Ues aux piopri6tes physiques dW milieu complexe tel que le sous-sol, 
rep6i6s k partir de donn6es obtenues par exploration du sous-sol ou par des mesures in situ, 
5 (6v6nements d'ordre g6ologique, g^ophysique, 6venements li6s k des donnees de 
production de fluides etc.), ces 6v6nements 6tant rep6r6s par des points dans un espace 
multivariable et d^finis par un vecteur d'attributs. 

Un autre exemple est donn6 par le brevet FR 01/05.675 du demandeur, qui decrit 
une mSthode pour faciHter le suivi au couis du temps de revolution d'une zone souterraine 
10 par analyse comparee d'un certain nombre n de jeux d'enregistrements sismiques obtenus 
successivement h. Tissue de n campagnes sismiques 3D successives (sismique dite 4D), 
dans laquelle on utilise une technique de reconnaissance des formes appliqute k I'ensemble 
des 6venements sismiques de plusieurs campagnes, consid6r6s et analyses simultan6ment. 

Citons aussi le brevet EP -A-671 017 (US 5 638 269) du demandeur, qui d6crit une 
15 m6thode permettant de rapprocher des doraifies g6ologiques obtenues par carottage ou 
diagraphies dans des puits et des donn6es sismiques obtenues par une campagne 
d'exploration sismique qui s'appuie sur une technique de calibrage statistique avec un 
etalonnage statistique obtenu en rapprochant les donn6es geologiques locales (mesur6es 
dans des puits) avec des attributs sismiques lus sur des traces sismiques obtenues au 
20 voisinage imm^diat de chaque puits. 

Les attributs sismiques qui font I'objet de cette interpretation sont calcul6s dans le 
domaine apres sommation, ou dans le domaine avant sommation. 

Dans le domaine apr^s sommation, les attributs classiques sont calculus i partir des 
amplimdes au niveau du reservoir, ou k partir de l'imp6dance P estim6e k I'aide d'une 
25 inversion strati graphique. 

Dans le domaine avant sommation. le nombre d'attributs pertinents pent 
consid6rablement augmenter. Dans ce cas, les attributs sont soit les amplitudes i diff6rents 
offsets (d6ports) ou angles d'incidence, soit des paramfetres issus d'une inversion 
stratigraphique jointe avant sommation (Imp6dances P et S, masse sp6cifique, produit des 
30 paramfetres de Lame avec la masse sp6cifique, etc.). 



Une methode dite d'analyse en composantes principales (ACP) bien connue des 
gens de I'art, peut etre appliqu6e k ces attributs afin d' analyser les relations existantes entre 
eux et de reduire le nombre d'attributs significatifs. Les composantes principales extraites 
de r ACP sent les nouveaux attributs : Us definissent une base vectorielle orthogonale ou 
5 non (rotation) et sont des combinaisons lindaires de variables initiales. lis peuvent atre 
utiUs6es dans des algorithmes statistiques de reconnaissance de foime et coirespondent k 
un filtrage multivarie de rinfonnation sismique initiale. Di£f6centes ^plications de 1' ACP 
sont d€crites> par exemple, dans : 

Dumay, J., Foumier, P., 1988, "Multivariate statistical analyses applied to seismic facies 
10 recognition". Geophysics, 53, 1 151-1 159 ; 

Hagen, D., C, 1982, "The application of principal components analysis to seismic data 
sets" : Geoexpl., 20, 93-1 11 .; ou encore 

dans le brevet FR-A- 2 772 138 (US 6 345 108) du demandeur. 

Mais I'ACP ne prend pas en compte la notion de groupe d'attributs et, par 
15 consequent, a est souvent difficile de donner une interpretation physique aux dbmposantes 
principales, ou bien de les reUer clairement aux attributs initiaux, surtout si deux-ci sont 
organises en groupes de mSme signification physique. 

L'analyse dite canonique est aussi une m6thode statistique bien connue des gens de 
I'art, qui permet d'^tablir les relations pouvant exister entre deux groupes de variables afin 
20 de savoir si ces deux groupes d&rivent les mSmes propri^tes. Cette methode est ddcrite par 
exemple par : 

- Hotelling, H., 1936, "Relations between two sets of variables", Biometrika, 28, 321- 
377. 

Un exemple d'application de l'analyse canonique est d^crite par exemple par : 

25 - Foumier. F., and Derain, J.F., 1995, "A statistical methodology for deriving reservoir 
properties from seismic data". Geophysics, 60, 1437-1450, 



Cette methode se limite k 1' etude de deux groupes de variables et done ne 
s'applique pas de fagon g6nerale a rinformation sismique multicubes. De plus, elle d6finit 
des variables synthetiques dans chacun des sous-espaces associ^s aux deux groupes, et non 
une seule base vectorielle permettant de compacter I'ensemble des variables initiales et de 
5 d6ciire les relations entre les deux groupes. 

Plusieurs m6thodes permettant de g^neraliser T analyse dite canonique en 6tendant 
I'analyse k plus de deux groupes de variables existent. Differents aspects de I'art anterieur 
dans le domaine consid6t6 sont d^crits, par exeraple, dans les publications suivantes : 

- Horst, P., 1961, "Relations among m sets of measures", Psychometdka, 26, n*2, 129- 
10 149 ; 

- Carol!, J.D., 1968, "A generalisation of canonical correlation analysis to three or more 
sets of variables", Proc. 76* Conv. Amer. Psych. Ass ; 

- Kettenring, J.R., 1971, "Canonical analysis of several sets of variables", 58, 3, 433- 
450 ;ou 

15 - Saporta, G., 1990, "Probabilit^s, analyse des donn^es et statistiques" , Technip, Paris, 

Ces m6thodes pr6sentent certaines limites. En particulier, elles ne permettent pas de 
d6crire chaque groupe s6par6ment, les variables synth6tiques ne peuvent 6tre relives 
facilement aux diff6rents groupes, et , par consequent, Stre physiquement interpr6t6es. De 
plus, seule la part d'inertie globale qu'elles repr6sentent peut etre connue, et non la part 
20 d'inertie qu'elles reprgsentent pour chaque groupe. E n'est par consequent pas possible de 
filtrer chaque groupe independamment. 

Une technique d'analyse 6galement connue dite analyse en composantes principales 
gen6ralis6e (ACPG), permet de comparer diff6rents groupes de variables tout en d6crivant 
chacun d'eux. On la trouve mise en oeuvre sur des donndes k caractere economique par 
25 exemple dans : 

- Casin, Ph, 2001, "A generalization of principal component analysis to K sets of 
variables", Computational Statistics & Data Analysis, 35, 417-428. 



La m^thode sdon-Pinventioii 

La m^thode selon I'invention propose une approche oil en exploitant une technique 
statistique multivariable appliquee & plusieurs jeux de traces d'enregisttement sismiques 
3D disponibles, obtenus par exploration sismique d'une zone souterraine, on compacte et 
on fatre les 6v6nements sismiques lisibles sur ces jeux de traces. La methode comporte : 

- la repartition de ces 6venements sismiques dans uii certain nombre p de famiUes de 
variables sismiques ayant chacune une signification physique particuli6re ; 

- la formation par combinaison des dites variables sismiques, de variables synthdtiques 
en nombre trfe inf^rieur au nombre des dites variables, que Ton obtient par 
construction d'une base vectorieUe orthogonale dans chacun des p ensembles d'analyse 
constitu6s par les donndes de chacune des dites families, d'oil d6coule la formation 
d'une base vectorieUe orthonorm^ d&rivant les p ensembles d'analyse ; 

- I'utiUsation de cette base vectorieUe orthonoimfie pour filtrer et d€criie les dits 
6v6iements sismiques. 

La formation des variables synthdtiques est obtenue par exemple en combinant une 
technique d'analyse en composantes principales ou ses techniques 66hv6es bashes sur la 
diagonaUsation d'une norme (recherche de vecteurs propres) et une technique de 
regression simple. 

Suivant un mode de mise en ceuvre, on peut soumettie les variables synth6tiques ^ 
une rotation pour am61iorer les corrSlations entre variables initiales et variables 
synth^tiques. 

Les variables synthetiques peuvent Stre form^es k partir de families comportant par 
exemple un ou plusieurs cubes de traces sismiques obtenues k Tissue d'une campagne 
sismique 3D, chacun correspondant k un mfime deport ou offset. 

On peut 6galement fonner des variables synth6tiques k partir de famiUes 
constitutes de cubes de traces sismiques obtenus k Tissue d'une campagne sismique 3D. 
chacun correspondant k un meme angle d'incidence. 



On peut encore foimer des variables synthetiques ^ partir de families constitutes de 
cubes d'attributs sisnjiques foim6s par traitement des amplitudes de traces sismiques 
obtenues ^ Tissue d'une campagne sismique 3D. 

■ On peut encore former des variables synth6tiques h partir de families constitutes de 
cubes de traces sismiques obtenues k Tissue de campagnes sismiques successives de la 



zone. 



On peut encore former des variables synthttiques a partir de families constitutes de 
cubes d'attributs sismiques formes par traitement des amplitudes de traces sismiques 
obtenues k Tissue de campagnes sismiques successives de la zone. 

10 La mtthode est particulidrement avantageuse en ce qu'elle peimet, appliqute k de 

I'infonnation sismique : 

- de dtfinir de nouveaux attributs sismiques, combinaisons lintaires des attributs 
sismiques initiaux, rtsumant Tinformation sismique initiale et permettant de filtrer 
cette information, tout en lespectant la structure multicubes des donntes 

15 - de montrer si tous les cubes d»attributs sismiques initiaux, ou seulement certains 
d'entre eux, apportent une information redondante, et de mesurer la part de cette 
information redondante 

- d'apporter une inteiprttation physique aux nouveaux attributs sismiques. 

Les nouveaux attributs sismiques ainsi formes peuvenl ensuite Stre analyses afin 
20 d'extraire les proprittts gtologiques de Tinformation sismique. Deux types de mtthode 
peuvent classiquement gtre utilists k cette fin : une premiere qui utiUse des techniques 
d'estimation statistique et foumit une interpittation quantitative en teimes de proprietts 
gtologiques, et une deuxi&me qui ufiUse des techniques de reconnaissance de forme pour 
analyser les faci&s sismiques et foumir une interpretation qualitatitive. 

25 Des exemples de ces deux types de mtthode sont mis en ceuvre respectivement 

dan^ les brevets EP-A-671 017 (US 5 638 269) et FR-A-2 768 818 (US 6051651) du 
demandeur, d6jk cit6s. 



Presentation succincte des figures 

les figures lA ^IC montrent I'organisation des enregistrements ' sismiques en p 
groupes : chaque groupe rassemble dans un cube I'ensemble des portions de traces 
conespondant k un critere physique identique ; 

les figure 2A k 2C montrent respectivement les cartes de trois premiers nouveaux 
attributs sismiques Z^^\ et ZP^ ; 

les figures 3A k 3D montrent respectivement le premier plan temporel vrai (3A) de la 
fengtre d*analyse et le mSme plan reconstitug/filtnS avec 1 (3B), 3 (3C)ou 8 (3D) des 
nouveaux attributs sismiques pour le cube d'imp6dance P ; 

les figures 4A k 4D montrent respectivement le premier plan temporel vrai (4A) de la 
fenatre d'analyse et le mSme plan reconstitu6/filti€ avec 1 (4B), 3 (4C)ou 8 (4D) des 
nouveaux attributs sismiques pour le cube d'imp6dance S ; 

les figures 5A k 5D montrent respectivement le premier plan temporel vrai (5A) et le 
mgme plan reconstitud/filtr6 avec 1 (5B), 3. (5Qou 8 (5D) des nouveaux attributs 
sismiques pour le cube de masse sp€cifique /? ; 

les figures 6 montrent respectivement les sections verticales, passant par le puits W2, 
vraie (6A), reconstito6e/filtree avec les 8 premiers nouveaux attributs (6B), ainsi que 
les diffdTMices (rfisidus) (6C) entre ces deux sections pour le cube d'impMance P ; 

-les figures 7 montrent respectivement les sections verticales, passant par le puits W2, 
vraie (7 A), reconstitu6e/filtr6e avec les 8 premiers nouveaux attributs (7B), ainsi que 
les differences (rfisidus) (7C) entre ces deux sections pour le cube d'imp6dance S ; 

la figure 8 montre la localisation des traces d'apprentissage utilis6es pour calibrer la 
fonction de classement dans le cadre de la reconnaissance des formes supervisee des 8 
premiers nouveaux attributs sismiques ; 

les figures 9A, 9B montrent respectivement les attributions des evenements sismiques 
et les attributions les plus fiables de ces mSmes 6venements dans le cadre de la 
reconnaissance des formes supervisee des 8 premiers nouveaux attributs sismiques ; et 



- les figures lOA, lOB represeatent respectivement les attributions des evenements 
sismiques et les attributions de ces memes evenements sismiques dans le cadre de la 
reconnaissance des formes non supervis6e des 8 premiers nouveaux attributs sismiques. 

Description detaill^e 

5 On dispose de plusieurs jeux d'enregistrements sismiques obtenus k des instants 

diff6rents par exploration sismique d^une zone souterraine etudiee. Les 6v6nements 
sismiques que I'on cherche k compactei/filtrer sont caract6ris6s^par des attributs sismiques 
souvent organises en groupes ou famiUes k signification physique particulifere (figures 1 A I 
IC). Ces attributs peuvent etre la succession des valeurs des amplitudes de portions de 

10 traces sismiques ou de paramfetres 61astiques issus d'une inversion 61astique jointe avant 
sommation, obtenue par exemple k I'aide de la methode decrite dans le brevet FR 2 800 
473 (EP 1 096 270)d6jk cit6 le long de la fenStre temporelle 6tudiee. 

1 - Notations 

Soit X une matrice compos6e de p tableaux juxtapos6s Xi,. . . Xj . . . Xp : 

X ={Xj ... X^ ... Xp}, X^=[Xp ... Xj.^ ... Xj,^,] 
15 oii Xj est un tableau compos6 de mj variables Xjj. de n individus (au sens statistique du 
tenne). La matrice X rassemble tous les 6v6nements sismiques disponibles. Ces. 
6v6nements sont tout d'abord organises en p groupes ou families Xj (figure lA). Chaque 
famiUe Xj correspond h une propriety physique particuli&re. Elle se pr6sente sous la forme 
d'un cube de valeurs coirespondant k cette prcpri6te, rep6rees par leur position spatiale 
20 Oigne. cdp, temps/profondeur). (figures IB et IC). Ces cubfes peuvent 6tre form6s par 
exemple par un ensemble des valeurs d'amplitude de n traces ou portions de traces 
assodtes k un deport commun ou St un angle commun ou encore la distribution spatiale des 
valeurs d'un attribut commun d6iv6 de ces amplittides telle que, par exemple, TimpSdance 
acoustique en ondes P ou S, la masse sp6cifique ou des param^tres de Lam6, dans une 
25 fenStre temporelle choisie de taille mj. Dans cet exemple, chaque variable X^.ic rassemble 
les valeurs de I'attribut sismique consid^re (amplitude sismique ou propri6t6 61astique) 
pour un ^chantiUon temporel k I'interieur de la fenStre d'analyse. Dans cet exemple, un 
cube comprend autant d' individus que de traces disponibles. 



Les mj variables Xj,k d'une table Xj definissent un espace denote Sj et sont 
supposfies centimes. Soit Dn la matrice diagonale des poids des individus. 

Soit un op6rateur de projection orthogonale sur T espace Sj defini par : 

nj^Xj(x'jD„XjYx'jD„ 

2 - Calcul des nouveaux attributs 

Le but de TAnalyse en Composantes Principales G6n6ralis6e (ACPG) est k la fois 
de decrire les proximit^s entre les p families de variables et de d6crire s6parement chaque 
groupe grace h des variables synth6tiques qui sent des combinaisons Uneaires des variables 
initiales. Alors, ce qui est coimnun k tons les groupes ou families, ou a certains d' entre eux 
seulement, et par cons6quent ce qui est different, pent etre mis en evidence. L*ACPG 
permet de r6duiie le nombre de variables utilise pour d6crire un jeu de donn6es et permet 
un filtrage multivarifi tout en conservant la structure muMcubes des donn6es initiales. Afih 
d'atteindie ce but, la premiere variable synth6tique calcul6e, que Von nomme Z^^\ calcul6p 
doit Stre telle que : 

!• les p projections orthogonales Z^^^ de Z^^^ sur les espaces Sj doivent Stre le plus proche 
les unes des autres ; 

2. Dans chaque groupe j la dispersion des n individus projetes sur Z^^^ est maximum. 

1. max y\R\Z^'\Zf') 

2. max |;Cov^(ZS^X,J 

Ceci se traduit mathematiquement par les conditions suivantes : 

ou est le coefficient de correlation multiple et Zj^^^ est la projection de Zf^^ sur Tespace 
Si norm^e k 1. 
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Ensuite,- on cherche une nouvelle variable synth6tique Z^^^ selon les memes criteres. 
On ajoute une condition d'orthogonalite entre les projections les variables Z^^ du mSme 
groupe, Le processus s'airete quand une base orthogonale dans chaque espace Sj est 
trouv6e. Le nombre de variables Z^^ n'excfede pas la dimension la plus grande des espaces 
Sj. 

Des methodes de rotation bien connues des gens de Tart, peuvent etre noises en 
oeuvre dans chacun des espaces ainsi trouv&, afin d'optimiser les coefficients de 
correlations entre les nouvelles variables et les variables initiales. 

La premifere variable Z^^^ est vecteur propre norme de la matrice de covariance 
XX^Dn correspondant k la plus grande valeur propre. Z^^^ est ensuite projete sur les p 
espaces Sj et, dans chaque espace, la regression des mj variables par Zj^^^ est calcul6e. La 
matrice des p rdsidus juxtaposes X/^^ d6finit le nouvel espace d'6tude. Ce nouvel espace 
correspond k Tespace initial duquel Tinfluence de Zj^^^ a et^ retir6 dans chaque groupe Sj. 
Les mj variables de la table Xj^^ d6finissent un sous espace Sj^^ de Sj. La seconde 
composante est le vecteur propre nonn6 de la matrice de covariance X^^^X^^^^ des 
r^sidus correspondant ^ la valeur propre la plus grande, tP^ est projetee sur les p espaces 
Sj^^\ Et la regression des mj variables par Zj^^^ est calculee dans chaque espace Sj^^\ La 
regression effectuee dans chaque espace Sj assure que la projection de Z^^^ sur Tespace Sj^^^ 
est orthogonale k la projection de Z^^^ sur T espace Sj, 

Le calcul des vecteurs propres est effectue par toute methode connue des gens de 

Tart. 

En pratique, la matrice XX'Dn etant d'ordie n, n pouvant etre tres grand, k une etape 
k, on calcule le vecteur propre norme V^^ de la matrice X^^^D^X^^ associe k la plus grande 
valeur propre. Le vecteur Z^^^ est ensuite simplement calcule par transformation lineaire 
Z<^>=X^>V^>et norme. 

Chaque groupe peut 8tre reconstitue avec un petit nombre d*attributs. Les variables 
initiales du groupe j peuvent s'exprimer comme la sonrnie d'une combinaison lineaire des 
variables Zj^^ et d'un residu. Si toutes les variables Zj^^ d'un groupe sont utilisees, le 
residu est nul. Les coefficients de la combinaison lineaire ne sont autres que les 
coefficients calcuies k Tissu de chaque regression dans les espaces Sj. 
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La m6thode selon invention a pour propriety, entre autres, que les variables Z^^ sont 
orthogonales deux k deux. 

En vue de TintCTpretation des variables synthetiques calculees, on calcule ensuite 
pour chaque groupe j : * ' 

5 1- le carr6 de la correlation entre la variable Z^^ et la variable Zj^^ du groupe j, indiquant 
le degr6 de proximit6 entre ces deux variables; 

2. le pourcentage d'inertie du groupe j expliqu6 par la variable Zj^\ 

Les con:61ations entre les variables Z^^ et les variables initiales Xjjc sont dgalement 
calcul6es pour mettre en 6vidence d'6ventuelles relations entre les Xj,k . Les correlations 
10 entre les variables 2^^^ et Xj,k d*im meme groupe j, ou les correlations entre variables Zj^^ 
pour un mSme k peuvent egalement Stre calculees. 

Les nouvelles variables synth6tiques ainsi cr66es constituent les nouveaux attributs 
sismiques recherches. 

3 - Example d'application ^ 

15 La m^thode selon invention a 6x6 appliqu6e h. trois cubes de propri6tes 61astiques 

(impedance P not6e Ip, impedance S not6e Is et masse sp6cifique notee p) obtenus aprds 
inversion elastique jointe avant sommation k Taide de la m6thode decrite dans le brevet 
ER2 800 473 d6ja cite d'enregistrements sisnadques regroup6s en cubes k angles 
dMncidence commuh. Chaque propriety elastique engendre un groupe de variables ou 

20 d* attributs sismiques. La fenStre temporfeUe d' analyse se compose de 19 echantillons de 
signaux pr61ev6s sur les enregistrements sismiques. Nous consid6rons comme variables les 
valeurs de la propri6t6 61astique pour un 6chantillon temporel k TintSrieur de la fenStre 
d' analyse. Chaque groupe est alors compos6 de 19 variables compos6es d'autant 
dMndividus que de traces disponibles. Comme il y a trois groupes de variables, nous 

25 disposons au total de 57 variables. Les trois cubes de paramfetres ont 6t6 prealablement 
horizontalises par rapport k une surface de reference. 

Les variables synthetiques Z^^ obtenues r6sument les informations respectives des 
trois propriet6s 61astiques, tout en prenant en compte les relations entre elles. 



ler aepoi 



3.1 - Analyse des relations entre les impedances S et la masse sp^cifique p 

Le tableau 1 ci-dessous domie la proportion de variance en pourcentage expliquee 
par les projections des 8 premieres variables synth6tiques Z^'^ k Z^^^ dans les trois 
groupes. n montre qu'il suffit de 8 variables pour rfisumer 87,4% de la variance de 
rimp6dance P, 94,4% de la variance de rimp6dance S et 88,4% de la variance de la masse 
sp6cifique. 







2(2) 


2^(3) 


2;(4) 










• 


Ip 


23,4 


24,5 


6.3 


13,1 


9,8 


3.1 


3,8 


3,4 


87,4 


Is 


26,6 


14,2 


15,2 


12,6 


5,7 


7,2 


5,8 


7,1 


94.4 


• 


16,9 


14.1 


13,6 


10,2 


8,0 


9,4 


10,1 


6,1 


88,4 



Le tableau 2 ci-apr6s donne le carr€ de la correlation entre les variables 
synthetiques et leurs projections sur les trois groupes. Cette valeur indique le degr6 de 
proximite entre la variable et le groupe. Plus la valeur est proche de 1, plus la variable est 
proche de sa projection et done du groupe. 











2(4) 


2;(5) 








Ip 


0,48 


0,79 


0,38 


0.72 


0,71 


0,24 


0,02 


0.30 


Is 


0,91 


0.43 


0,86 


0,53 


0,36 


0,59 


0,09 


0,73 


• 


0,78 


0,49 


0,77 


0,46 


0,65 


0,79 


0,96 


0,46 
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Par exemple, Z^^^ est tres proche de Is et proche de p , mais tres 61oigne de Ip. 
Alors que Zf^^ est proche de Ip, mais trds 61oign6 de Is et de p. Sur les S premieres 
variables trois sont proches de Ip seul, deux sont proches k la fois de Is et de>o , une 
est proche de Is seul, et deux sont proches de p seul. Ainsi, il apparalt que Ip et Is ne 
5. contiennent pas la iiiSme information et qu*une part de Tinfonnation contenue dans Is et 
dans p semble commune. 

3.2 - Filtrage des impedances P, S et de la masse specifique 

Les trois propriety 61astiques peuvent etre reconstitutes avec un nombre restreint 
de nouveaux attributs sisnaiques- Le nombre de variables synthetiques utilis6es pour la 
10 reconstitution pent etre different d*un groupe k Tautre : il d6pend du pourcentage d'inertie 
que Ton desire reconstituer/filtrer. 

3.3 - Recomiaissance des formes supervisee et non supervisee des nouveaux 
attributs r 

n est possible d'interprtter geologiquement les nouveaux attributs sismiques au 
15 travers d'une analyse de facifes avec un algorithme de reconnaissance de forme supervis6 
ou non. 

En reconnaissance des formes supervisee, il est n&essaire de choisir des classes 
d'apprentissage regroupant des 6v6nements sismiques situ6s aux voisinages de puits dont 
r interpretation g6ologique a €l€ e£fectu6e. Cette m6thodoIogie permet de classer 
20 Tensemble des 6v6nements sismiques en fonction des caractSres g6o]ogiques observes aux 
puits. 

Dans notre cas, les connaissances g6ologiques sur les trois puits disponibles Wl, 
W2 et W3, ainsi que les connaissances geologiques g6n6rales sur la zone d' etude nous 
amenent k considerer trois classes : une premi&re classe correspondant aux puits Wl et W2 

25 representent des sables turbiditiques, une seconde classe correspondant au puits W3 
represente des boues turbiditiques et une troisieme classe represente des h6mip61agites. Les 
puits Wl et W2 montrant des caract6ristiques geologiques proches, seuls les ev6nements 
sismiques situes autour du puits W2 sont utilises, le puits Wl servant de test aveugle. La 
localisation du puits W3 6tant incertaine, les 6venements sismiques utilises pour 

30 caracteriser la classe 2 sont situ6s autour d'un puits virtuel d6fini ^ proxinciite de W3. 
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Enfin, ne disposant d'aucun puits pour la classe 3, un puits viituel est defini k partir de nos 
connaissances geologiques g6n6rales. La Fig.7 montre la localisation gfiographique des 
trois classes d'apprentissage. Chaque classe est constitu6e de 121 6v6neinents sismiques 
(11x11). Ces ev6nements servent k calibrer une fonction de classement qui est ensuite 
appliqude h Tensemble des evenements sismiques. 

Les Fig.8A, 8B montrent respectivement les affectations obtenues et les affectations 
les plus probables, c'est-i-dire ayant une probabilite d* affectation sup6rieure ou egale a 
80%. Le puits Wl, utilis6 en test aveugle, est correctement attribu6 k la meme classe que le 
puits W2. Le puits W3 est localise en bordure des classes 1 et 2. La classe 2 se distribue le 
long d'une fonne chenalisante* La carte de la figure 8B montre que les resultats sont 
stables : peu d*6v6nements presentent une mauvaise probabilite d' affectation. * 

En reconnaissance des formes non supervis6e, les classes d'apprentissage sont 
constitutes en utilisant les modes d*une fonction de density de probability multivarite 
calcul6e k partir de Tensemble des 6v6nements sismiques, dans un espace engendre par les 
huit premiers nouveaux attributs. Les dvenements sismiques de plus forte density de 
probabilite constituant les pics de cette fonction de densit6 sont s61ectionnes pour realiser 
Tapprentissage de la fonction de classement. Cette fonction est ensuite appliquee k 
Tensemble des 6v6nements sismiques consideres pour les affecter k Tune des classes : les 
Fig.9A, 9B montrent les affectations obtenues, ainsi que les affectations les plus probables, 
c'est-lL-dire ayant une probability d' affectation sup6rieure ou egale k 80%. Les cartes 
obtenues sont tn^s proches des cartes des fig. 8A, 8B, La classe 3 de T analyse supervisee 
correspond aux classes 3, 4 et 5 de Tanalyse non supervisee. Les puits Wl et W2 sont 
attribu6s k la mSme classe et le puits W3 est toujours situ6 k la fipontifere entre les classes 3, 
et surtout 1 et 2. Cette convergence des resultats valide rinterpr6tation de Tapproche 
supervisee, 
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REVENDICATIONS 

1) M6thode pour compacter et filtrer des evenements sismiques 3D lisibles sur 
plusieurs jeux de traces d'enregistrement sismiques disponibles, obtenus par exploration 
sismique d\ine zone souterraine, caracteris6e en ce qu'elle comporte : 

- la r6partition de ces 6v6nements sismiques dans un certain nombre p de families de 
variables sismiques ayant chacune une signification physique particuliere ; 

- la formation par combinaison des dites variables sismiques, de variables synth6tiques 
en nombre txes inferieur au nombre des dites variables, que Ton obtient par 
constmction d*une base vectorielle orthogonale dans chacun des p ensembles d'analyse 
constitu6es par les donn6es de chacune des dites fanailles, d'ou ddcoule la formation 
d'xme base vectorielle orthonorm6e decrivant les p ensembles d'analyse ; et 

- Tutilisation de cette base vectorieUe orthonorm^e pour filtrer et dScrire les^ dits 
6v6nements sismiques. 

'"V 

2) M6thode selon la revendication 1, caract6ris6 en ce que la formation, des 
variables synth6tiques.est obtenue en combinant une technique d'analyse en composantes 
principales ou des m6thodes d6riv6es de diagonalisation d'une norme et une technique de 
regression simple. 

3) Methode selon la revendication 1 et 2, caract6risee en ce que Ton somnet k une 
rotation les variables synth6tiques pour am^liorer les correlations entre variables initiales et 
variables synth6tiques. 

4) Methode selon la revendication 1 ou 2 ou 3, caracteris6e en ce que Ton forme des 
variables synth6tiques k partir de families constituees de cubes de traces sismiques obtenus 
a Tissue d'une campagne sismique 3D, chacun correspondant k un meme deport ou offset 

5) M6thode selon la revendication 1 ou 2 ou 3, caract6risee en ce que Ton forme des 
variables synth6tiques a partir de families constitu6es de cubes de traces sismiques obtenus 
k Tissue d'une campagne sismique 3D, chacun correspondant k un mSme angle d'incidence. 
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6) Methode selon la revendication 1 ou 2 ou 3, caract6risee en ce que Ton forme des 
variables synthetiqpies a partir de families constituees de cubes d'attributs sismiques 
obtenus par traitement des amplitudes de traces sismiques obtenues k Tissue d'une 



7) M6thode selon Tune des revendications prec6dentes, caracterisee en ce que Ton 
forme des variables synth6tiques a partir de families constituees de cubes de traces 
sismiques obtenus par des explorations sismiques successives de la zone. 

8) M6thode selon Tune des revendications precedentes, caracterisee en ce que Ton 
forme des variables synth6tiques a partir de families constituees de cubes d'attributs 
sismiques form6s par traitement des amplitudes de .traces sismiques obtenues par des 
explorations sismiques successives de la zone. 
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9 Class 3 : hemipelagjte 
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